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Abstract, The exact treatment of nuclei starting from the constituent nucleons and the fundamental interactions 
among them has been a long-standing goal in nuclear physics. Above all nuclear scattering and reactions, which 
require the solution of the many-body quantum-mechanical problem in the continuum, represent an extraordinary 
theoretical as well as computational challenge for ab initio approaches. We present a new ab initio many-body ap- 
proach which derives from the combination of the ab initio no-core shell model with the resonating-group method 
[4]. By complementing a microscopic cluster technique with the use of realistic interactions, and a microscopic 
and consistent description of the nucleon clusters, this approach is capable of describing simultaneously both 
bound and scattering states in light nuclei. We will discuss applications to neutron and proton scattering on s- 
and light p-shell nuclei using realistic nucleon-nucleon potentials, and outline the progress toward the treatment 
of more complex reactions. 
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1 Introduction 

To understand the evolution of the Universe, we need to 
understand nuclear reactions. Indeed, low-energy fusion 
reactions represent the primary energy-generation mech- 
anism in stars, help determining the course of stellar evo- 
lution, and are crucial in the formation of the chemical ele- 
ments. In addition, much of what we know about neutrino 
oscillations is established from neutrino emerging from the 
Sun following the j6-decay of reactions products, particu- 
larly '^B. The light-ion fusion reactions which encompass 
the Standard solar model need to be understood better, if 
solar neutrinos are to provide even more precise informa- 
tion on the neutrino oscillation properties. As an example, 
the ^Be(/?, yf'^ radiative capture is a rather poorly-known 
step in the nucleosynthetic chain leading to '^B, which in 
turn is the dominant source of the high-energy solar neu- 
trinos (through jS-decay to '^Be) detected in terrestrial ex- 
periments. 

Furthermore, nuclear reactions are one of the best tools 
for studying exotic nuclei, which have become the focus of 
the next generation experiments with rare-isotope beams. 
These are nuclei for which most low-lying states are un- 
bound, so that a rigorous analysis requires scattering con- 
ditions. In addition, much of the information we have on 
the structure of such short-lived nuclei is inferred from re- 
actions with other nuclei. 

Unfortunately, the calculation of nuclear reactions rep- 
resent also a formidable challenge for nuclear theory, the 
main obstacle being the treatment of the scattering states. 
In this paper we will present a brief overview of existing 
theoretical methods for nuclear reactions, highlighting in 
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particular recent progress in the ab initio calculation of 
low-energy scattering of light nuclei. 



1 .1 Overview of reaction approaches 

Because of their importance nuclear reactions attract much 
attention, and there have been many interesting new devel- 
opments in the recent past. In this section we will give a 
brief overview of the theoretical efforts devoted to nuclear 
reactions, and in particular scattering of light nuclei. How- 
ever, this is not intended to be completely exhaustive. 

Nuclear reaction approaches may be classified accord- 
ing to two broad categories. The first category embraces 
the so-called microscopic approaches, in which all the nu- 
cleons involved in the scattering process are active degrees 
of freedom, and the antisymmetrization of the many-body 
wave functions is treated exactly. 

In the three- and four-nucleon sectors there has been 
remarkable progress in the past ten years: the Faddeev [11, 
Faddeev-Yacubovski |2 3 1, Alt-Grassberger and Sandhas 
(AGS) 14151, hyperspherical harmonics |6|, Lorentz inte- 
gral transform methods 1718191 . etc., are among the best 
known of several numerically exact techniques able to de- 
scribe reactions observables starting from realistic nucleon- 
nucleon (A^A^) and in some cases also three-nucleon (NNN) 
forces. 

Going beyond four nucleons there are fewer ab initio 
or ab initio inspired methods able to describe reactions 
observables starting from realistic forces. Only very re- 
cently the Green's function Monte Carlo ifTOl . the no-core- 
shell model combined with the resonating group method 
(NCSM/RGM) 0111121 and the fermionic molecular dy- 
namics lfT3l have made steps in this direction. 
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Reactions among light nuclei are more widely described 
starting from semirealistic A^A^ interactions with adjusted 
parameters within the traditional resonating-group method 
II14I15I16I17I18I or the generator coordinate method II19I20I21I . 
which are microscopic cluster techniques. 

A second category is that of few-body methods de- 
scribing scattering among structureless clusters. Here one 
starts form nucleon-Nucleus (usually optical) potentials fit- 
ted on some reaction observable, and the nucleus core is 
usually inert. There are exact techniques which can em- 
ploy either local or non-local optical potentials like the 
Faddeev or AGS methods adopted by Deltuva II22I23I . and 
various approximated ones, hke the continuum-discretized 
coupled channel equations II24I25L distorted wave born ap- 
proximations, or various adiabatic approximations 1261 . etc., 
which usually adopt local optical potentials. 

Finally, there are also some recent attempts to describe 
reactions among light nuclei in an effective-field theory ap- 
proach for halo nuclei |27 28|. Starting from experimental 
resonance parameters for the system under investigation, 
phase shifts and cross section are predicted at low energy. 

Here we focus on one of the above mentioned approaches, 
i.e. the NSCM/RGM, for which we will present formalism 
in Sec. |2] and a collection of results in Sec. [3] Conclu- 
sions and an outline of possible future developments will 
be given in Sec.|4] 



2 Formalism 

A brief overview of the NCSM approach is presented in 
Sec. 12. II the resonating-group method is introduced in Sec. 12.2 
and the NCSM/RGM formalism is described in Sectionl23] 



2.1 Ab initio no-core shell model 



techniques based on the second quantization and devel- 
oped for standard shell model calculations can be utilized. 
As a downside, one has to face the consequences of the 
incorrect asymptotic behavior of the HO basis. 

Standard, accurate A^A' potentials, such as the Argonne 
V18 (AV18) 1 29 1, CD-Bonn |30|, INOY (inside non-local 
outside Yukawa) [31 1 and, to some extent, also the chiral 
N^LO I.32J . generate strong short-range correlations that 
cannot be accommodated even in a reasonably large HO 
basis. In order to account for these short-range correla- 
tions and to accelerate convergence with respect to the in- 
creasing model space, the NCSM makes use of an eff'ec- 
tive interaction obtained from the original, realistic A^A^ or 
NN+NNN potentials by means of a unitary transformation 
in a n-body cluster approximation, where n is typically 2 
or 3 |33|. The effective interaction depends on the basis 
truncation and by construction converges to the original 
realistic A^A^ or A^A^ + NNN interaction as the size of the 
basis approaches infinity. 

On the other hand, a new class of soft potentials has 
been recently developed, mostly by means of unitary trans- 
formations of the standard accurate A^A^ potentials men- 
tioned above. These include the Viowk 1 34], the Similarity 
Renormalization Group (SRG) |35| and the UCOM |36| 
A^A^ potentials. A different class of soft phenomenological 
A^A^ potential used in some NCSM calculations are the so- 
called JISP potentials |37|. These so-called soft potentials 
are to some extent already renormalized for the purpose 
of simplifying many-body calculations. Therefore, one can 
perform convergent NCSM calculations with these poten- 
tials unmodified, or "bare." In fact, the chiral N^LO A^A^ 
potential |32| can also be used bare with some success. 
NCSM calculations with bare potentials are variational with 
respect to the HO frequency and the basis truncation pa- 
rameter A^max- 



The NCSM is a technique for the solution of the A-nucleon 
bound-state problem. All A (point-like) nucleons are active 
degrees of freedom, hence the difference with respect to 
standard shell model calculations with inert core. Starting 
from a microscopic Hamiltonian (p,- being the momentum 
of the ith nucleon and m the nucleon mass) 
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2.2 Resonating-group method 



The resonating-group method (RGM) II14I15I16I17I18I38I 
is a microscopic cluster technique in which the many-body 
Hilbert space is spanned by wave functions describing a 
system of two or more clusters in relative motion. Here, we 
will limit our discussion to the two-cluster RGM, which 
is based on binary-cluster channel states of total angular 
momentum J, parity n, and isospin T, 



containing realistic A^A^ (Vfi'^) or A^A^ plus A^A^A^ (Vfij^'^) 
forces (both coordinate- and momentum-space interactions 
can be equally handled), the non-relativistic Schrodinger 
equation is solved by expanding the wave functions in terms 
of a complete set of A-nucleon harmonic oscillator (HO) 
basis states up to a maximum excitation N,„a.xfiQ above the 
minimum energy configuration, with Q the HO frequency. 
The choice of the HO basis in a complete Nmaxfi^^ space 
is motivated by its versatility. Indeed, this is the only basis 
which allows to work within either Jacobi relative coordi- 
nates or Cartesian single-particle coordinates (as well as 
easily switch between the two), while preserving the trans- 
lational invariance of the system. Consequently, powerful 
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In the above expression, {\A-aaiI'^^Tij and |aa2^^'7"2) 
are the internal (antisymmetric) wave functions of the first 
and second clusters, containing A-a and a nucleons (fl<A), 
respectively. They are characterized by angular momentum 
quantum numbers /i and h coupled together to form chan- 
nel spin s. For their parity, isospin and additional quantum 
numbers we use, respectively, the notations jij, T,-, and a,, 
with i - 1,2. The cluster centers of mass are separated by 
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the relative coordinate 



, A-a , A 

TA-fl.a - rA-aM^A-cua - -. / , l"; / , I" ; , (3) 

A- a ^ a . ■f-' , 

;=1 j=A-a+l 

where {r,, / = 1, 2, ■ ■ ■ , A) are the A single-particle coordi- 
nates. The channel states (|2|i have relative angular momen- 
tum {. It is convenient to group all relevant quantum num- 
bers into a cumulative index v = [A-a a\I"^ T\ ; a 02!^ T2', 
si}. 

The former basis states can be used to expand the many- 
body wave function according to 
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However, to preserve the Pauli principle one has to intro- 
duce the appropriate inter-cluster antisymmetrizer, schemat- 
ically 



^v — 



(A-fl)!fl! 
A! 



J](-rp, 



(5) 



where the sum runs over all possible permutations P that 
can be carried out among nucleons pertaining to differ- 
ent clusters, and p is the number of interchanges char- 
acterizing them. Indeed, the basis states ^ are not anti- 
symmetric under exchange of nucleons belonging to dif- 
ferent clusters. 

The coefficients of the expansion (|4]i are the relative- 
motion wave functions gl'^ir). These are the only unknowns 
of the problem, to be determined solving the non-local 
integral-differential coupled-channel equations 



Yj^drr" [^i:jir', r) - E Ni:j ir' , r)] 
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= 0, 



(6) 
where E denotes the total energy in the center-of-mass 
frame. Here, the two integration kernels, specifically the 
Hamiltonian kernel, 

H'7(r', r) = (<;| X'Hm. \<Pi7) , (7) 

and the norm kernel. 



Nt:jir', r) = (0^:;| ^,-j^v K^) , 



(8) 



contain all the nuclear structure and antisymmetrization 
properties of the problem. The somewhat unusual presence 
of a norm kernel is the result of the non-orthogonality of 
the basis states (|2|, caused by the presence of the inter- 
cluster antisymmetrizer. The exchange terms of this anti- 
symmetrization operator are also responsible for the non- 
locality of the two kernels. 

The main inputs of the RGM method are i) the internu- 
cleon interaction; and //) the wave functions of the (A - a)- 
and a-nucleon clusters. Staring from the latters, in the tra- 
ditional RGM the clusters internal wave functions are of- 
ten (but not exclusively) translationally invariant HO shell- 
model functions of the lowest configuration or a linear su- 
perposition of such functions. The value of the HO size 



parameters b are chosen ad hoc to reproduce properties of 
the nucleon clusters (such as size and/or binding energy, 
etc.) within the adopted interaction. This somewhat sim- 
plified description of the clusters internal wave functions 
is usually compensated by the use of semirealistic A^A^ in- 
teractions, such as the VoUcov ||39l or Minnesota BOll po- 
tentials, with parameters that can be adjusted to reproduces 
important properties of the compound nucleus or reaction 
under study. The spin-orbit force, not present in the men- 
tioned semi-realistic interactions, is sometimes added to 
the microscopic Hamiltonian. Exception to this general de- 
scription of the RGM approach exist, particularly in the 
few-nucleon sector, where the method has been utilized in 
combination with realistic A^A^ and NN+NNN forces |38|. 
Finally, the treatment of the Coulomb interaction between 
charged clusters does not represent an issue in the RGM 
approach. 

The advantage in expressing the RGM basis states (|2]) 
as antisymmetrized products of single-particle functions, 
and in particular Slater determinants, lies in the ability to 
carry out analytical derivations of the required matrix ele- 
ments (|7]i and ^. Once the integration kernels are calcu- 
lated, by solving the integral-differential coupled channel 
equations (|6|l subject to appropriate boundary conditions, 
one obtains bound-state wave functions and binding en- 
ergies or scattering wave functions and scattering matrix, 
from which any other scattering and reaction observable 
can be calculated. 



2.3 Ab initio NCSM/RGM approach 



A new first-principles, many-body approach capable of si- 
multaneously describing both bound and scattering states 
in light nuclei has been developed by combining the RGM 
with the ab initio NCSM 0111121 . This new approach com- 
plements the microscopic-cluster technique of the RGM 
with the utilization of realistic interactions and a consistent 
microscopic description of the nucleonic clusters, while 
preserving important symmetries such as Pauli exclusion 
principle, translational invariance, and angular momentum. 
More in detail, the formalism presented in Section l272l can 
be combined with the ab initio NCSM as follows. 

First, we note that the Hamiltonian can be written as 



i/=r,,i(r)+n^,ei + Vc(r) + i/, 
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H(a 
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where H(^A-a) and //(„) are the {A— a)- and a-nucleon intrin- 
sic Hamiltonians, respectively, TieKr) is the relative kinetic 
energy and 'Viei is the sum of all interactions between nu- 
cleons belonging to different clusters after subtraction of 
the average Coulomb interaction between them, explicitly 
singled out in the term Vc{r) - ZxyZ^vC^ jr (Ziy and Z2v 
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being the charge numbers of the clusters in channel v) 
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Nuclear, V^Cr—rj, cr,, ctj, t,, tj), and point-Coulomb com- 
ponents of the two-body potential have been listed explic- 
itly ( (Ti, Ti denoting spin and isospin coordinates, respec- 
tively, of the ith nucleon). If the A-nucleon Hamiltonian 
contains a A^A^A^ force, the inter-cluster interaction T^rei 
will present also a contribution from the latter, denoted 
here with Vf/^^ ,. 

The cluster's Hamiltonians and inter-cluster interaction 
'Vrei are consistent, as they contain the same realistic po- 
tentials. The clusters internal wave functions are also treated 
consistently: ( \A-a a\I"^ Tn and \a ffi^^'T^i) are obtained 
by diagonalizing H^A-a) and H^a), respectively, in the model 
spaces spanned by the NCSM basis. Both the (A - a)- and 
fl-nucleon model spaces are characterized by the same HO 
frequency Q and maximum number A^max of excitations 
above the minimum configuration. At the same time, in 
calculating the Hamiltonian and norm kernels of Eqs. (|7]l, 
and ([8]), all "direct" terms arising from the identical permu- 
tations in both ^y and ^y> are treated exactly (with respect 
to the separation r) with the exception of / (P;^,'^ | 'Vi-ei \^i'/)- 
The latter and all remaining terms are localized and can be 
obtained by expanding the Dirac 6 of Eq. (|2l) on a set of 
HO radial wave functions with identical frequency Q, and 
model-space size A^max consistent with those used for the 
two clusters. In this respect, we note that, thanks to the 
subtraction of the average potential Vcir), 'Viei is localized 
also in the presence of the Coulomb force. 

If the adopted potential generates strong short-range 
correlations, we employ consistent NCSM effective inter- 
actions derived from it. More specifically, the cluster eigen- 
states are obtained by employing the usual NCSM effective 
interaction 1 33 1 . However, in place of the bare potential en- 
tering 'Vrei we adopt a modified effective interaction, which 
avoids renormalizations related to the kinetic energy. Fol- 
lowing the notation of Ref. |33|, at the two-body cluster 
level this is given by V^^g - Hi^s - ^2eff' where H'^^^ is 

//02 + Vi'2, 

Vn in the limit 



the effective Hamiltonian derived from H^ 

'^lA 



-mQ~i 



with V'.j — —mi,j-r-/A. inolb mai v, « — * Vff 
Nmax — > °° and, for each model space, the renormaliza- 
tions related to the kinetic energy and the HO potential in- 
troduced in H2eS are compensated by the subtraction of 
^2e«- ^^^ kinetic-energy renormalizations are appropri- 
ate within the standard NCSM, but they would compro- 
mise scattering results obtained within the NCSM/RGM 
approach, in which the relative kinetic energy and average 



Coulomb potential between the clusters are treated exactly 
(that is, are not truncated within a finite HO model space). 



2.3.1 Integration kernels 



To give a somewhat more in depth description of the for- 
malism involved in the calculation of the matrix elements 
O and ([8]i, here we will present examples of algebraic ex- 
pressions derived within the single-nucleon projectile ba- 
sis, i.e., for binary-cluster channel states (|2]i with a - \ 
(with channel index v = {A-1 ci/^Ti; 1 \\; st]). In this 
model space, the norm kernel is rather simple and is given 
by 



Ni:j{r',r) = {cPi:j,\ 1 -Y,p,A KO 



(11) 



= 6v'v 



6(r' - r) 



r' r 
{cPIZ\Pa-x.aH7) , 



-{A-l)Y^Rn,Ar')RnM 
(12) 



where it is easy to recognize a direct term, in which initial 
and final state are identical (corresponding to diagram (a) 
of Fig. [1), and a many-body correction due to the exchange 
part of the inter-cluster anti-symmetrizer (corresponding to 
diagram {b) of Fig. [U. As the exchange Pa-i,a is a short- 
range operator, in calculating its matrix elements we re- 
placed the delta function of Eq. (|2| with its representation 
in the HO model space. Such HO expansion is appropriate 
whenever the operator is short-to-medium range. 

The presence of the inter-cluster anti-symmetrizer af- 
fects also the Hamiltonian kernel, and in particular the ma- 
trix elements of the interaction: 



A-l 



U'7('-'' '■) = (^vv I -^f 1 - Z ^A I^v7) (13) 

1=1 

= f f rei(r') + Vc{r')+E'^, ' ^' ] Ni:j(r' , r) 

(14) 



+^°(r',r)+%%(r',r), 



where £ ', 'are the eigenenergies of the (A - l)-nucleon 

"1 
cluster. If no A'A'A^ forces are present in the Hamiltonian 

one obtains a "direct" term involving interaction and ex- 
change of one of the nucleons in the first cluster with the 
nucleon (a = 1) of the second cluster (see diagrams (c) and 
(d) of Fig.[TJ, and an "exchange" term involving the inter- 
action of the Ath nucleon with one of the (A - 1 ) nucleons, 
accompanied by the exchange with a second of such nucle- 
ons. Diagram (e) of Fig.[T]describes this latter term. These 
two potential kernels, which together constitute the matrix 
element ('^y"^\ "y^i^ |^vr^/' have the following expres- 
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1 2-A-2A 






(c) 



{d) 



(e) 



Fig. 1. Diagrammatic representation of: (a) "direct" and (b) "exchange" components of the norm kernel; (c and d) "direct" and (e) 
"exchange" components of the potential kernel for the (A - 1,1) cluster basis. The first group of circled lines represents the first cluster, 
the bound state of A-1 nucleons. The separate line represents the second cluster, in the specific case a single nucleon. Bottom and upper 
part of the diagram represent initial and final states, respectively. 



sions: 



n'n 
x{0iZ\VA-Ul-PA-lA)H7) 



(15) 



%%(/, r) = -(A - 1)(A - 2)Y,Rn'e'{r')R,Ar) 

n'n 

X (0;:;| Pa-i,a Va-i,a-x K7)- (16) 

The inclusion of a A^A^A^ interaction in the Hamiltonian 
is straightforward, and amounts to extra "direct" and "ex- 
change" potential kernels, which can be obtained in a sim- 
ilar way. 

Being translationally-invariant quantities, the Hamilto- 
nian and norm kernels (|7] [8]) can be "naturally" derived 
working within the NCSM Jacobi-coordinate basis. How- 
ever, particularly for the purpose of calculating reactions 
involving /?-shell nuclei, it is computationally advantageous 
to introduce Slater-determinant (SD) channel states of the 
type 



(sT) 



l'fv7>SD = [{\A-aa,hT,)5o\aa2hT2)) 

^ye{R'^^l)t''^^Rnt{RtL), (17) 

in which the eigenstates of the (A-a)-nucleon fragment 
are obtained in the SD basis (while the second cluster is 
still a NCSM Jacobi-coordinate eigenstate), and R*°' = 
^-1/2 2^^^^ J r, is the vector proportional to the center of 
mass coordinate of the a-nucleon cluster Indeed, it easy to 
demonstrate that translationally invariant matrix elements 
can be extracted from those calculated in the SD basis of 
Eq. ( fTTb by inverting the following expression: 



sd(<;|4..K7)sd 



'„n,e,,j. 



n',t'„n,l,,J, 
X 



x{n,.trNU\OOnee)^^ {nXNLe\00n'{'{')^ ■ (18) 



Here (3t.i. represents any scalar and parity-conserving and 
translationally-invariant operator (Ot.i = .^, ^H^, etc.), 
and {nr{rNL{\QOn{{)^, {n%NU\OQn'£'£'}^ are general 
HO brackets for two particles with mass ratio a/(A -a). We 
exploited both Jacobi-coordinate and SD channel states to 
verify our results. 

As an example, the single-nucleon projectile "exchange" 
part of the norm kernel within the Jacobi-coordinate basis 
for a system of A = 3 nucleons is given by: 
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Hem{nieiSiIiTi\2aiI1'Ti)iind{n[{[s[I[T[\2a\I^'T[)iim 
the coefficients of the expansion of initial and final two- 
nucleon target wave functions, respectively, with respect 
to the HO basis states depending on the Jacobi, spin, and 
isospin coordinates ^i = (r i - ri)/ "Vz, cri , 0-2, and ti , T2, 
respectively. 



{^io-iO-2TiT2\nieiSiIiTi} , 



(21) 



where n\,{i are the HO quantum numbers corresponding 
to the harmonic oscillator associated with ^1, while si,Ii, 
and Ti are the spin, total angular momentum, and isospin 
of the two-nucleon channel formed by nucleons 1 and 2, re- 
spectively. Note that the basis (1211 1 is anti-symmetric with 
respect to the exchange of the two nucleons, (-)'^i+'i+^i = 
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r' [fm] 



5^ 



Fig. 2. Calculated "exchange" part of the norm kernel, 
Nff^ir' , r) - [,]. for the n + a(g.s.) 'S 1/2 (top panel), and ^P3/2 
(bottom panel) channels as a function of the relative coordinates 
r and r', using the N^'LO A'A' potential 1 32] at tiQ= 19 MeV. The 
"5 1/2 channel is strongly influenced by the Pauli-exclusion prin- 
ciple, which forbids to accommodate more than four nucleons 
into the i-shell of a nuclear system. The four nucleons forming 
the "'He g.s. sit mostly in the OHQ shell. Accordingly, in the ^S 1/2 
channel the "exchange"-part of the norm kernel suppresses the 

(dominant) OhQ contribution to the d function of Eq. J12t (and, 

1+ 1 
consequently, to the 5 -wave relative-motion wave function gj^^^ ) 

coming from the fifth nucleon in i-shell configuration. The situ- 
ation is difi'erent if we consider a different partial wave, like the 
^/'3/2 channel, corresponding to the 'He resonance. In this case 
the exchange norm represents a small and positive correction to 
the delta function. 



of the exchange operator Pa-\,a yields: 
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„TfO^'''^\\A-lahTi)sD are 
one-body density matrix elements of the target nucleus and 
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n'f'j'-m' ^-iii'j 



Next we extract 



the corresponding translationally-invariant matrix elements, 
i-i\<l>^y^n!'^^''''^}, by inverting Eq. dH for 
Pa-i,a- The final step follows easily from 



^ V,. n,. ' 

a = 1 and Ot 
Eq. (O. 

Due to the exchange terms of the intercluster antisym- 
metrizers, norm and potential kernels are non-local and ap- 
pear as surfaces in three dimensions such as, e.g., those 
shown in Figs. |2] and [3] The latter figures present results 
of single-channel calculations carried out using n-a cluster 
channels with the a particle in its g.s. (note that the index 
y = {4 g.s. 0^0; Ij j', jC] is simply replaced by the quan- 
tum number £). The interaction models adopted are the 
N^LO A^A^ potential |32|, and the Vio^k NN potential [341 
derived from AV18 1291 with cutoff A = 2.1 fm"'. 



2.3.2 Orthogonalization of the RGM equations 

An important point to notice, is that Eq. (jSJ does not repre- 
sent a system of multichannel Schrodinger equations, and 
gl'^ir) do not represent Schrodinger wave functions. This 
feature, which is indicated by the presence of the norm 
kernel N:^'Jir',r) and is caused by the short-range non- 
orthogonality induced by the non-identical permutations in 
the inter-cluster anti-symmetrizers, can be removed by in- 
troducing normalized Schrodinger wave functions 



x; 



'^(r) 



-^f 



dy y Nyy(r, y) 



y 



(23) 



-1. Finally, {n'e',n[£[,A\nifun{,Ah are the general HO 
brackets for two particles with mass ratio 3. 



At the same time, as mentioned above, the matrix el- 
ements of the operators Pa-i,a, Va-i,a(1 - Pa-\,a)' and 
PA-i.A^A-2,A-i '^an be more intuitively derived working 
within the SD basis of Eq. SvH . Using the second-quantization 
formalism, they can be related to linear combinations of 
matrix elements of creation and annihilation operators be- 
tween (A-l)-nucleons SD states. As an example, the case 



where N^ is the square root of the norm kernel, and ap- 
plying the inverse-square root of the norm kernel, N^^, 
to both left and right-hand side of the square brackets in 
Eq. (|6]l. By means of this procedure, known as orthogo- 
nalization and explained in more detail in Ref. |12|, one 
obtains a system of multichannel Schrodinger equations: 

Itdr) + Vc(r) -(E- eH'J' - fif ^^)]^^ 



■J]Jdr'r'^w!7(r, 



,.x7(r') 



r) 



= 0, 



(24) 
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Fig. 3. Calculated "direct" (top figures) and "exchange" (bottom figures) potential kernels for the n + Q'(g.s.) 5 1/2 (left column), and 
^Pi/i (right column) channels as a function of the relative coordinates r and r', using the Via„k NN potential 134] at TiQ = 18 MeV. The 
^Si/2 channel is strongly influenced by the Pauli-exclusion principle, which forbids to accommodate more than four nucleons into the 
i-shell of a nuclear system: the exchange potential introduces repulsion between a nucleon in the a particle and the fifth nucleon when 
they are both in s-shell, largely suppressing the attractive contribution of the direct potential. The situation is different if we consider 
a different partial wave, like the ^P3/2 channel, corresponding to the 'He resonance. In this case both direct and exchange potential are 
attractive. 



f'T- 

where E^^. are the energy eigenvalues of the /-th clus- 
ter (/ - 1,2), and W^"J(r', r) are the overall non-local po- 
tentials between the two clusters, which depend upon the 
channel of relative motion, while do not depend upon the 
energy E of the system. 



3 Results 

The two-cluster NCSM/RGM formalism within the single- 
nucleon projectile outlined in the previous section, can be 
used to calculate nucleon-nucleus phase shifts below three- 
body break threshold, by solving the system of multi- 
channel Schrodinger equations ( l24l i with scattering bound- 
ary conditions. In the next sections we will review part 
of the results for neutrons scattering on ^H, ^He and '"Be 
and protons scattering on ^■''He, using realistic A^A^ poten- 
tials, which were first presented in Refs. 1. 11 J and I.12J . and 
present some new calculations. 



3.1 Convergence with respect to the HO model 
space 

To study the behavior of our approach with respect to the 
HO model space, we have performed NCSM/RGM scatter- 
ing calculations for the A - 5 systems, using the Viowt NN 
potential [341, which is "soft" and we treated as "bare", 
and the N^LO A^A^ interaction [32 1, which generates strong 
short-range correlations, thus requiring the use of effective 
interactions. In particular, for this convergence tests, we re- 
stricted our binary-cluster basis to target-nucleon channel 
states with the target in its g.s. (corresponding to channel 
indexes of the type v = {4 g.s. 0+0; l^^^; \t}). 

Results obtained for Viow* are presented in Fig.|4] The 
overall convergence is quite satisfactory, with a weak de- 
pendence on A^max- 

Figure |5] presents the convergence rate (achieved by 
using two-body effective interactions tailored to the HO 
model-space truncation) obtained for the same n-a scat- 
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Fig. 5. Dependence on #„,„ of the n-a(g.s.) phase shifts with the 
N'LO 132J A'A' potential at fiO = 19 MeV. 



tering phase shifts with the N^LO potential. Clearly, the 
N-'LO results converge at a much slower rate than the Viow*: 
ones. However, a gradual suppression of the difference be- 
tween adjacent A^max values with increasing model-space 
size is visible, although the pattern is somewhat irregular 
for the P phase shifts. 

Although not shown, the p-a phase shifts present anal- 
ogous convergence properties. 

The next figure. Fig. |6] compares the N,mx = 17 re- 
sults for the previously discussed Viowi and N^LO A^A^ in- 
teractions, and those obtained with the CD-Bonn A^A^ po- 
tential |30|. The NCSM/RGM calculations for the latter 
potential were carried out using two-body effective inter- 




E,^ [MeV] 

Fig. 6. Calculated n-a{g.s.) phase shifts obtained in the largest 
model space (A'max = 17) for the Viow* I.34J A' A' potential at fiQ = 
18 MeV and the N^LO |32) and CD-Bonn 03 NN interactions 
at no =19 MeV. 



actions, and present a convergence pattern similar to the 
one observed for N^LO. Clearly, the ^Pi/2 and ^^3/2 phase 
shifts are sensitive to the interaction models, and, in par- 
ticular, to the strength of the spin-orbit force. This obser- 
vation is in agreement with what was found in the earlier 
study of Ref. I.IOJ . Following a behavior akeady observed 
in the structure of p-shell nuclei, CD-Bonn and N^LO in- 
teractions yield about the same spin-orbit splitting. On the 
contrary, the larger separation between the ViowA 3/2" and 
1/2" resonant phase shifts is direct evidence for a stronger 
spin-orbit interaction. 

As the 1 /2^ channel is dominated by the repulsion be- 
tween the neutron and the a particle induced by the Pauli 
exclusion principle (see also Figs. |2] O, the short-range 
details of the nuclear interaction play a minor role on the 
^S 1/2 phase shifts. As a consequence, we find very similar 
results for all of the three adopted A^A^ potential models. 



3.2 Test of the NCSM/RGM approach in the 
four-nucleon sector 

A stringent test-ground to study the performances of the 
NCSM/RGM approach within the single-nucleon projec- 
tile basis is provided by the four-nucleon system. Numer- 
ically exact calculations for the A - A sector have been 
already successfully performed within accurate few-body 
techniques, such as Faddeev-Yacubovsky, AGS, and hyper- 
spherical harmonics methods II2I3I4I5I6I41I . 

Figures|7]and|8]show n -^Y\ and p -■'He phase shifts, re- 
spectively, calculated with the N^'LO A^A^ potential |32 1 to- 
gether with results obtained by Deltuva and Fonseca |4 43 
from the solution of the AGS equations (+ symbols) us- 
ing the same interaction. The convergence behavior of The 
NCSM/RGM calculations was achieved using two-body 
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of the H nucleus (i.e. the channels v = {3 g.s. 



\; si]). 



effective interactions tailored to the model-space trunca- 
tion, as outlined in Sec. 12.31 For the '^o, '^ i and ^S i par- 
tial waves, the increase in model-space size produces grad- 
ually smaller deviations with a clear convergence towards 
the A^max = 19 results. The rest of the phase shifts, particu- 
larly the ^Pq, show a more irregular pattern. Nevertheless, 
in the whole energy-range we find less than 2 deg absolute 
diff'erence between the phases obtained in the largest and 
next-to-largest model spaces. 

Concerning the comparison to the highly accurate AGS 
results, in general the agreement between the two calcu- 
lations worsens as the relative kinetic energy in the cm. 
frame, £],_[„, increases. This discrepancy is a manifestation 
of the the influence played by closed channels not included 
in our basis states, that is, target-nucleon channel states 
with the target above the /[' = | g.s., and 2h-2 config- 
urations, both of which are taken into account by the AGS 
results. Because these states coiTespond to the breakup of 
the A = 3 system, it is not feasible to include them in the 
current version of the NCSM/RGM approach, which so far 
has been derived only in the single-nucleon projectile ba- 
sis. However, we are planning on extending our approach 
to be able to account for the target breakup, and these de- 
velopment will be discussed in future publications. Never- 
flieless, this test of the NCSM/RGM approach in the A = 4 



sector clearly indicates that one has to pay attention not 
only to the convergence with respect to the HO model- 
space size Nm-dx, but also to the convergence in the RGM 
model space, which is enlarged by including excited states 
of the nucleon clusters in the binary-channel basis states. 



3.3 ''He(Af,Af)''He scattering 

A better scenario for the appHcation of the NCSM/RGM 
approach within the single-nucleon projectile basis is the 
scattering of nucleons on "^He. This process is character- 
ized by a single open channel up to the "^He breakup thresh- 
old, which is fairly high in energy. In addition the low- 
lying resonances of "^He are naiTow enough that they can 
be reasonably reproduced by diagonalizing the four-body 
Hamiltonian in the NCSM model space, and consistently 
included as closed channels in the NCSM/RGM model 
space. 

In Fig. |9] we explore the effect of the inclusion of the 
first six excited states of the "^He on the n-a scattering 
phase shifts obtained with the N^LO A^A^ interaction. More 
specifically, in addition to the single-channel results (dot- 
ted line) discussed in Sec. 13.11 we show coupled-channel 
calculations for five different combinations of '^He states. 
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nucleus (i.e. the channels v = {3 g.s. 



^^1). 



i.e., g.s.,0+0 (dash-dotted line), //) g.s.,0+0,0 (dash- 
dot-dotted line), //Og.s.,0+0,0"0, 1^0, 11 (dash-dash-dotted 
line), iv) g.s.,0+0, 2^0 (dashed line), and u) g.s.,0+0, 2^0, 2^1 
(solid line). 

The use of these five different combinations of ground 
and excited states (also shown in the legends of Fig. ^ 
indicates that the ^5 1/2 phase shifts are well described al- 
ready by coupled channel calculations with g.s. and first 
0^0 (the ^S 1/2 phase shifts obtained in the four larger Hilbert 
spaces are omitted for clarity of the figure). On the con- 
trary, the negative parity excited states have relatively large 
effects on the P phase shifts, and in particular the 0, 1 
and 1^1 mostly on the ^Pi/2, whereas the 2^0 and 2^1 
on the ^P3/2- These negative parity states influence the P 
phase shifts because they introduce couplings to the i-wave 
of relative motion. Though also /^' = 1" couples to ^ = 
in the 3/2" channel, the coupling of the /^' = 2" states is 
dominant for the ^^3/2 phase shifts. 

Figure [10] compares an accurate /^-matrix analysis of 
the nucleon-Qf scattering |43 1 with NCSM/RGM results ob- 
tained including the first six '^He excited states as in Fig.|9] 
This comparison reveals that for both neutron (top panel) 
and proton (bottom panel) projectiles we can describe quite 
well the ■^S 1/2 and qualitatively also the ^^3/2 phase shifts, 
using the N^LO A^A^ potential. On the contrary, the same 
interaction is not able to reproduce well the two P phase 
shifts, which are both too small and too close to each other. 
This lack of spin-orbit splitting between the ^Pi/2 and ^P]/2 
results can be explained by the omission in our treatment 
of the AAA terms of the chiral interaction, which would 
provide additional spin-orbit force. This sensitivity of the 
P phases to the strength of the spin-orbit force corrobo- 



rated by the differences among the Viowt, N^LO and CD- 
Bonn results in Fig.H] ^Pi/2 and ^^3/2 are both larger and 
more apart for Viow* than for N^LO or CD-Bonn potentials. 
At the same time, the good agreement of the N^'LO and 
CD-Bonn ^S 1/2 phase shifts with their Viow/t analogous and 
the /^-matrix analysis owes to the repulsive action (in this 
channel) of the Pauli exclusion principle for short nucleon- 
a distances, which has the effect to mask the short-range 
details of the nuclear interaction. 

As shown in Fig. [TT] a stronger spin-orbit interaction 
is shown also by the '*He(n, n)*He 3/2" and 1/2" resonant 
phase shifts obtained using the SRG-evolved N^LO AA 
potential with cutoff A = 2.02 fm"' [44], in a NCSM/ 
RGM n-a model space including g.s. and 0^0 states of "^He. 
As for Viow* the two resonances present larger separation. 
In addition, in the SRG-N^LO case, the ^Pi/2 phase shifts 
lie on top of the n-a R-mathx analysis of Ref. |43 1 in the 
whole energy range, and one finds a satisfactory agreement 
with the /^-matrix results for all four lowest partial waves 
starting from about /iyn =12 MeV. This particularly good 
quality of the SRG-N^LO A-'^He phase shifts for energies 
far from the low-lying resonances is the reason behind the 
the fairly good agreement, presented in Fig.[T2] of the cal- 
culated *He(n, «)'*He analyzing power (top panel) and dif- 
ferential cross section (bottom panel) with the experimen- 
tal data of the 17 MeV polarized -neutron experiments of 
Ref. ||45i. 



3.4 n-i**Be scattering and "Be bound states 

Figure [13] highlights one of the promising aspects of the 
NCSM/RGM approach, that is the abihty (through the use 
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Fig. 9. Influence of the lowest six excited states (/[Ti = 
0+0,0^0, rO, ri, 2^0,2-1) of the a particle on the n-a ^Si,2 
(top panel), ^Pi/2 (central panel), and ~P3/2 (bottom panel) phase- 
shift results for the N^LO NN potential ^ at hO = 19 MeV. 
Dotted (g.s.) and dash-dotted (g.s., 0*0) lines correspond to 
single- and coupled-channel calculations in a A'.^ax = 17 model 
space, respectively. The effects on the ^Pi/2 and ^^3/2 phase shifts 
of the further inclusion of, respectively, the 0, 1^0, 1^1, and 
2^0, 2^1 states are investigated in a N^-jj. = 15 model space. 



of SD channel states) to perform ab initio scattering calcu- 
lations for p-shell nuclei. The ^S 1/2 (top panel) and ^^5/2 
(bottom panel) n-^^Be phase shifts were obtained in a Nmax - 
6, ti£2 =13 HO model space. The inclusion of the 2| ex- 
cited state of '"Be has a significant effect on the S and 
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spectively, of Fig. |9] The ^03/2 phase shifts were obtained in a 
coupled-channel calculation including ground and first 0* excited 
state of ■'He, in a Ny„.^^ = 17 HO model space. 



more importantly on the D phase, where it is essential for 
the appearance of a resonance below 3 MeV. We note that 
a resonance has been observed at ~ 1.8 MeV with a ten- 
tative spin assignment of (5/2, 3/2)* B6l . The further ad- 
dition of the 2^ and especially 1| excited states produces 
rather weak differences. We have also extracted the scat- 
tering length for the ■^S 1/2 partial wave, and found a result 
of -i-10.7fm, which is comparable to the value of -1-13.6 fm 
obtained by Descouvemont in Ref. |47|, by fitting the ex- 
perimental binding energy of "Be. 

Although we mainly described its scattering applica- 
tions, the NCSM/RGM is a powerful tool also for struc- 
ture calculations, particularly for loosely-bound systems. 
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Fig. 11. Calculated n-a phase shifts obtained in a N,^a^ = 17 
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N^LO (Id NN potential at W = 20 MeV, compared to an R- 
matrix analysis of data (+) 1431 . 



By imposing bound-state boundary conditions to the set 
of coupled channel Schrodinger equations of Eq. jTM . we 
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states of '"Be. The obtained ^S 1/2 scattering length is +10.7 fm. 



tested the performance of our single-nucleon projectile 
NCSM/RGM formalism for the description of one-nucleon 
halo systems. In particular, because of the well-known parity- 
inversion between its two bound states with respect to the 
predictions of the simple shell model |48|, the "Be nu- 
cleus represents an excellent test ground for our approach. 

Large-scale ab initio NCSM calculations with several 
accurate A^A^ potentials of the "Be low-lying spectrum were 
not able to explain its g.s. parity inversion |50|. The ex- 
planation for these results can be searched into two main 
causes: (i) the size of the HO basis was not large enough 
to reproduce the correct asymptotic of the «-'"Be compo- 
nent of the 11 -body wave function; (ii) the A^A^A^ force, not 
included in the calculation, plays an important role in the 
inversion mechanism. The second hypothesis was corrob- 
orated by the results obtained with the INOY A^A' poten- 
tial |31 1, which seemed to indicate the possibility to reach 
the inversion in a large NCSM basis. 

By studying the "Be bound states in a NCSM/RGM 
model space spanned by n-'^Be channel states with in- 
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Table 1. Calculated energies (in MeV) of the '"Be g.s. and of the lowest negative- and positive-parity states in "Be, obtained using the 
CD-Bonn NN potential [301 at ^Q = 13 MeV. The NCSM/RGM results were obtained using n+'^Be configurations with A'^a^ = 6 g.s., 
2|,2+,and 1| states of '"Be. 





A/ 

^ * max 


'"Be 


"Be(i 


) 


"Be(j") 




£g.,s. 


E 


E,i, 


E E,i, 


NCSM 149 501 
NCSM 149 50 511 
NCSM/RGM 151 1 
Expt. 


8/9 

6/7 


-57.06 
-57.17 

-64.98 


-56.95 
-57.51 
-57.59 
-65.16 


0.11 
-0.34 
-0.42 
-0.18 


-54.26 2.80 
-54.39 2.78 
-57.85 -0.68 
-65.48 -0.50 



elusion of the N„,ax - 6 g.s., 2|,22, and 1| of '"Be, we 
are now in the position to address the first hypothesis. In- 
deed, the coiTect asymptotic behavior of the n-'^Be wave 
functions is described naturally in the NCSM/RGM ap- 
proach. The energies of the lowest 1 /2^ and 1 /2" states 
of "Be obtained in the NCSM and in the NCSM/RGM 
calculations, using the same CD-Bonn NN interaction 1301 
at hQ =13 MeV adopted in Ref. Il50ll . are presented in Ta- 
ble [T] The relatively small differences between the A^max = 
6/7 and A^max - 8/9 NCSM results, seems to indicate a 
reasonable degree of convergence for these calculations. 
The 1/2" state appears to be the g.s., and the 1/2^ state 
is about 2.8 MeV above the n -"^Be threshold. A compari- 
son to the NCSM/RGM calculations (obtained in a model 
space including g.s., 2|, 22, and 1| states of "^Be) shows 
a rough agreement for the 1/2" state, whereas for the 1/2"'' 
state one observes a dramatic difference (~3.5 MeV) in the 
energy. The 1/2" and 1/2+ NCSM/RGM states are both 
bound and the 1/2+ state is the g.s. of "Be. Correspond- 
ingly, we obtain a B(E1; ^ ~* 5 ) value of 0.18 e^ fm^, 
which is not far from experiment. 

To understand the mechanism which makes the 1/2+ 
state bound in the NCSM/RGM, we evaluated mean val- 
ues of the relative kinetic and potential energies as well 
as the mean value of the "^Be energy, and compared them 
to those obtained by restricting all the integration kernels 
to the HO model space (i.e. by replacing the delta func- 
tion of Eq. (fT2l i with its representation in the HO model 
space). These results are shown in Table |2l The model- 
space-restricted calculation is then similar, although not 
identical, to the standard NCSM calculation. In particular, 
as in the NCSM one loses the correct asymptotic behav- 
ior of the n-'"Be wave function. We observe that in the 
full NCSM/RGM calculation both relative kinetic and po- 
tential energies are smaller in absolute value. This is an 
effect of the re-scaling of the relative wave function in the 
internal region when the Whittaker tail is recovered. The 



Table 2. Mean values of the relative kinetic and potential energy 
and of the internal '"Be energy in the "Be 1/2+ ground state. All 
energies in MeV. NCSM/RGM calculation as in Table[T] See the 
text for further details. 



NCSM/RGM 


(Trd) 


(W) 


£['"Be(g.s.,ex.)] 


^tot 


Model Space 
Full 


16.65 
6.56 


-15.02 
-7.39 


-56.66 
-57.02 


-55.03 
-57.85 



difference is significantly more substantial for the relative 
kinetic energy than for the potential energy. As a result one 
obtains a dramatic decrease of the energy of the 1 /2+ state, 
which makes it bound and even leads to a g.s. parity in- 
version. This study shows that a proper treatment of the 
coupling to the n -'"Be continuum is essential in explain- 
ing the g.s. parity inversion. However, we cannot exclude 
that the A^A^A^ force plays a role in the inversion mecha- 
nism, not until accurate calculations with both A^A^A^ force 
and full treatment of the «-'"Be tail will be performed. 



4 Conclusions and Outlook 

We have reviewed the NCSM/RGM, a new ab initio many- 
body approach capable of describing simultaneously both 
bound and scattering states in light nuclei, by combining 
the RGM with the use of realistic interactions, and a mi- 
croscopic and consistent description of the nucleon clus- 
ters, achieved via the ab initio NCSM. In particular, we 
have outlined the formalism on which the NCSM/RGM is 
based, and given examples of the algebraic expressions for 
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Fig. 14. Calculated «-'-C phase shifts obtained within the 
NCSM/RGM in a A'max = 16 model space including g.s. and 
first excited states of '^C, using the SRG-N^LO NN potential 
of Ref (44) with cutoff" A = 2.66 fm"' and HO frequency 
fiQ = 24 MeV. The '^C wave functions were obtained within 
the importance-truncated NCSM I53I54I . 
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(a) 



(b) 



(c) 



Fig. 15. Diagrammatic representation of: (a) "direct", (b) "single-exciiange", and (d) "double-excliange" components of tlie norm kernel 
for tiie (A - 2, 2) cluster basis. The first group of circled lines represents the first cluster, the bound state of A-2 nucleons. The second 
couple of circled lines represents the second cluster, in the specific case a deuteron nucleus. Bottom and upper part of the diagram 
represent initial and final states, respectively. See also caption of fig.[T] 



the integral kernels within the single-nucleon projectile ba- 
sis, working both with Jacobi, and SD single-particle co- 
ordinate bases. As the spurious cm. components present 
in the SD basis were removed exactly, in both frameworks 
the calculated integral kernels are translationally invariant, 
and lead to identical results. Several analytical as well as 
numerical tests were performed in order to verify the ap- 
proach, particularly by benchmarking independent Jacobi- 
coordinate and SD calculations for systems with up to 5 
nucleons. 

Among the applications, we reviewed results for neu- 
tron scattering on ^H, '^He and '"Be and proton scattering 
on ■'•''He, using realistic A^A^ potentials. Our A = 4 scatter- 
ing results were compared to earlier ab initio calculations 
performed in the framework of the AGS equations, and the 
convergence properties of the NCSM/RGM approach with 
respect to the adopted model spaces were discussed in de- 
tail. For the A = 5 system, we found that all adopted A^A^ 
potentials provide a fairly good (in some case excellent) 
description of the S -wave phase shifts. On the contrary, the 
P-wave phase shifts that we obtained with any of the real- 
istic A^A^ potentials present both insufficient magnitude and 
splitting with respect to the /^-matrix analysis of the data. 
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Fig. 16. Calculated "exchange" part of the norm kernel, 
Nfi^{r' , r) - ^^T-^ for the rf(g.s.)-l-Qr(g.s.) channels as a function 
of the relative coordinate r at r' = 1 fm, using the SRG-N^LO 
A'A' potential |44| at tiQ = 20 MeV. Effects of the Pauli-exclusion 
principle are visible partcularly in the 1+0 channel. 



Concerning the calculations performed with the N^LO A^A^ 
potential, it is anticipated that the inclusion of the A^A^A^ 
terms of the chiral interaction would lead to an enhanced 
spin-orbit splitting, and recover the predictions of the R- 
matrix analysis. We also showed that the SRG-N^LO inter- 
action presents a larger spin-orbit strength than the N^LO 
A^A^ potential itself, and allows a fairly good description 
of '^He(n, n)'^He angular cross section and analyzing power 
for energies far above the low-lying resonances. An impor- 
tant aspect of the NCSM/RGM approach is its suitability 
for the description of loosely-bound systems, such as the 
"Be nucleus. Although we cannot exclude that, e.g. the 
A^A^A^ force plays a role in the inversion mechanism, we 
have demonstrated that a proper treatment of the coupling 
to the n-'°Be continuum leads to a dramatic decrease of 
the energy of the i state, which makes it bound and even 
leads to a g.s. parity inversion. 



Since the publication of the first results firi21, the 
NSCM/RGM approach has been applied to the description 
of nucleon scattering on several p-shell nuclei, such as ^Li, 
'^C, and '^O (see, e.g., Fig.[T4li. Key to these calculations, 
which will be published in a forthcoming paper ll52l . are 
two factors: (/) the ability of the NCSM/RGM to take ad- 
vantage of the powerful second quantization techniques, 
while preserving the translational-invariance symmetry of 
the system; and (//) the use of large HO model spaces (large 
A^max values) for the expansion of the clusters internal wave 
functions, and hence of the short-range parts of the integra- 
tion kernels. The use of the so-called importance-truncated 
NCSM ll53l for the description of the clusters internal wave 
functions is essential in achieving the second of these two 
points. This method, first introduced by R. Roth and R 
Navratil in 2007 ll53l for the calculation of ground states, 
then further developed in Ref. Il54l . and now extended to 
the calculation of excited states 1521 . allows to reach large 
HO model spaces for/?-shell nuclei by means of an a-priori 
selection of the most important NCSM basis states. Using 
importance-truncated clusters wave functions within the 
NCSM/RGM formalism allows us to reach for scattering 
calculations on heavier nuclei the same level of conver- 
gence obtained here for the A - A and A - 5 systems. 

The NCSM/RGM formalism has been also extended to 
include two-nucleon (deuteron) projectiles, and the inclu- 
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sion of the three-nucleon (triton and -^He) and four-nucleon 
("^He) projectiles are planned ahead. 

Calculations of deuteron-nucleus scattering are under- 
way. As an example, Fig.[T5]shows the diagrammatic rep- 
resentation of the matrix elements involved in the calcula- 
tion of the norm kernel within the deuteron-projectile ba- 
sis, and the "exchange" component for such kernel for the 
d{g.s.)+a(g.s.) system is presented in Fig. [16] 

In addition, the coupling of the single-nucleon and two- 
nucleon projectile basis will allow the first ab initio calcu- 
lation of the ^H{d, n)'^He fusion. 

Further, it is possible and desirable to extend the binary- 
cluster {A-a, a) NCSM/RGM basis by the standard A- 
nucleon NCSM basis to unify the original ab initio NCSM 
and NCSM/RGM approaches. This will lead to a much 
faster convergence of the many-body calculations compared 
to the original approaches and, most importantly, to an 
optimal and balanced unified description of both bound 
and unbound states. Extensions of the approach to include 
three-body cluster channels are also among our future plans, 
and the feasibility of such a project is supported by recent 
developments on the treatment of both three-body bound 
and continuum states (see, e.g., Refs. II21I55I56I57I58II ). 
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